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Abstract 

4^ ■ 

^ ' A novel Information Theory based method for determining the density of 



states from prior information is presented. The energy dependence of the 



^ ' density of states is determined from the observed number of states per energy 



interval and model calculations suggest that the method is sufficiently reliable 

to calculate the thermal properties of nuclei over a reasonable temperature 

range. 
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The simplest expression for the nuclear level density has been obtained in the Fermi-gas 
model by Bethe and later modified by Bloch There are, however, a number 

of shortcomings in this approach. For example, the lack of coupling to the collective part 
of the nuclear spectrum leads to an energy-independent level density parameter. Recently 
there has been considerable theoretical activity in the determination of the nuclear many- 
body density of states, taking into account shell, pairing, and deformation effects [^H, 
finite size effects [^, and thermal and quantal effects |1T0| , |TT| as well as improvements in 



the determination of the spin cut-off factors [T^. Furthermore, the multiple inverse Laplace 



transform used to determine the nuclear density of states from the Grand-Canonical partition 
function of a Fermi gas appears to lead to certain inconsistencies in the folding of nuclear 



level densities [13|. In spite of these deficiencies Bloch's formula is widely used, particularly 



as a means of parameterizing the experimentally determined nuclear level densities 
Generally the level density parameter is taken to be constant. Should it be energy dependent, 
an a priori form for the energy dependence must be assumed. 

In the present work we will demonstrate a simple means of extracting the nuclear level 
density from the experimental information, based on Information Theory |TB|. We wish to 



determine the distribution of states or level density, p{E), for a nucleus such that 

roo 

/ piE)Q{Ei - E)Q{E - Ei+,)dE = ^ = l,...,iV (1) 
Jo 

where Ni, the number of states in the energy interval -Ej+i — Ei, is assumed known, and 
B (x) = 1 if X > and is zero otherwise. As a zeroth-order guess we assume that the 
distribution of states or level density is given by the Bloch formula 



2JaoE 



(2) 



where Oq is the level-density parameter, often expressed in terms of the nuclear mass number 
A by 

A 

ao = (3) 

Kq 



with, for example, ko = 8 MeV. The actual distribution p{E) may now be determined by 
minimizing the information entropy of p{E) relative to Pq{E) 



S{p,po)= fp{E) In ^dE 



(4) 



subject to the experimental information given in eq. (|ip. This yields 



p{E) = po(^)exp 



N 



Y,X,Q{E,-E)Q{E-Ei+i] 
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(5) 



where the Lagrange multipliers Aj for the intervals Ei < E < -Ej+i are determined from the 
experimental information given in eq. (P. 

Now suppose we have determined the set of Aj and we put 



A(^.) = Xi 



(6) 



where Ei is the midpoint of the energy interval Ei < E < -Ej+i. If the X{Ei) vary smoothly 
with energy, a least-squares fit to this quantity should enable us to extrapolate the p{E) 
beyond the region in which the experimental information is known. If this is the case the 
density of states may now given by 



piE) = po{E)e-'^^\ 



(7) 



or, from eqs (|) and 



PiE) 
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exp 
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Note here that the energy dependence of the density of states is determined from the prior 
experimental data. 

Thermal properties of nuclei are calculated from the nuclear partition function Z (T) 
for fixed nuclear mass number A, where T is the nuclear temperature. The contribution to 
Z (T) of the continuous part of the nuclear spectrum is 



roo 

Z{T)= p (E) e-^'^dE 

■^-^min 



(9) 



where i?min is some suitable lower limit. 

To test the proposed method we have performed the following model calculations. Rather 
than use actual experimental information for the number of states Ni in eq. (|1|), we have 
generated the Ni by assuming that for energy E > 1 MeV the experimental distribution of 
states is given by 



2. 



AE 



(10) 



kiE)_ 

where k = k{E) to allow the energy dependence of the density of states to deviate from the 
Bloch form. The size of the energy intervals -Ej+i — Ei in eq. (0) was fixed at 0.4 MeV, 
with ten such intervals covering the energy range 1-5 MeV. For a particular choice k{E) the 
Lagrange multipliers Aj were calculated from eqs (0), (I), (D and (D, with A arbitrarily set 
to 56 in eq. (^. A series of least-squares polynomial fits of low order was then performed 
for the modified Lagrange multipliers 



X'(e,) = ^^. (11) 
^ ^ 2VAE ^ ' 

The many-body density of states and nuclear partition function were then calculated from 
eqs (H) and (H) , respectively. 

A sample fit to A' [Eij is shown in figure 1 for the choice 

k{E) = 8 (12) 

in eq. ([10|), corresponding to the Bloch formula (0) with /cq = 8; a perfect fit should produce 
A' = identically. Since we take the number of states to be an integer, whereas the integral 
in (|lD with the density of states given by ( p!OD produces a real number, we do not expect 
to find all the Lagrange multipliers to be zero. The error bars included in the figure were 



generated by arbitrarily changing the number of states to ± y Ni/2. They serve both 
to indicate the sensitivity of the calculation to uncertainties in the input, and also to provide 
a scale for A'. The fit was also weighted according to these error bars. As figure 1 indicates, 
the quality of the fit increases as the order of the polynomial is increased, but A' = 4 x 10~^ 
(constant) is already satisfactory. 



The corresponding many-body density of states was calculated from eq. with E 
extrapolated beyond the region of fit 1-5 MeV to 20 MeV; the density relative to the input 
density eq. (|T0[) is shown in figure 2. It is seen that A' = constant provides an almost perfect 
fit for all energies. Polynomial fits of higher order, whilst giving a slightly better fit up to 
E = 5 MeV, cannot be extrapolated much beyond this energy. 

The nuclear partition function was computed using eq. (H), with -Emin = 1 MeV. Due 
to the rapid increase in the calculated density of states for higher order polynomial fits to 
the modified Lagrange multipliers, the infinite upper limit of the energy integration was 
replaced by -Emax, determined as follows. First, ZiT) for a given temperature was calculated 
using the input density of states eq. (0)) with the upper limit gradually increased until 
the integral converged to a constant value. The final value for the upper limit, -Emaxj which 
varied between 30 MeV and 150 MeV depending on T, was then used in the calculation for 
that temperature with the fitted density of states eq. (H). The results are plotted in figure 

3, which shows Z{T) calculated with eq. relative to the partition function calculated 
with the input density of states. The fit A' = constant produces an almost perfect result up 
to T = 2 MeV, whereas the higher order fits are far from satisfactory except at the lowest 
temperatures. 

These calculations were repeated with other choices for k{E) in eq. 
the use of 

k{E) = 6 + E/2 (13) 

is illustrated in figures 4-6. This choice allows a reasonable energy dependence for the level 
density parameter, permitting the density of states to deviate sufficiently from the Bloch 
form to provide a realistic test of the method. (It should be noted that in practice the energy 
dependence of k{E) is constrained; too strong a dependence causes the density of states to 
increase rapidly at higher energies or to decrease suddenly to zero.) As seen from figure 

4, the lowest order polynomial which fits A' (^Eij within the error bars is linear in energy. 
Figure 5 shows that this same fit provides a satisfactory description of the density of states 
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10). For example. 



well beyond the energy region in which the fit was produced. The quality of the calculated 
partition function, shown if figure 6, is excellent. 

Finally, we have checked the calculations for sensitivity to various details in the method. 
For example we have performed calculations in which: 

• the energy Ei for the energy interval E'j+i — Ei was computed as the mean energy 
weighted according to the density of states, rather than simply as the midpoint energy; 

• the fit to A' (^Ei^ was not weighted according to the uncertainty in A'; 

• the fit was to A (^Ei^ rather than to A' (^Ej^ ; 

• the energy range 1-5 MeV was spanned by 5 intervals of width 0.8 MeV, or 16 intervals 
of width 0.25 MeV, rather than 10 intervals of width 0.4 MeV. 

In most cases the results were not changed significantly. There is a slight preference for 
fewer, wider energy intervals, and for the fit to A' to be weighted. A fit to A', rather than A, 
is definitely preferred, at least for those functions k{E) actually used. 

A simple Information theory based method for determining the density of states from 
prior experimental measurements has been presented. A good approximation of the density 
of states is obtained even at higher energies and any energy dependence is determined from 
the prior experimental data. However, the extrapolated density of states will not necessarily 
be very accurate if at higher energies processes occur which are not reflected in the lower 
energy prior information. With this proviso, the model calculations show that the partition 
function can be accurately calculated over a reasonable temperature range which suggest 
that the thermal properties of nuclei can be adequately described by this method. 
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FIGURES 

FIG. 1. Modified Lagrange multipliers for the case k{E) = 8. The circles indicate values 
computed with eq. (8 ), the solid line represents the fit with a polynomial of order zero, the dashed 
line the fit of order one and the dash-dot line the fit of order two. The calculation of the error bars 
is explained in the text. 

FIG. 2. Calculated density of states, relative to the exact density of states, for the case k{E) 
= 8. The solid line is computed with A fitted with a polynomial of order zero, the dashed line with 
the fit of order one and the dash-dot line with the fit of order two. 

FIG. 3. Calculated partition function, relative to the exact partition function, for the case k{E) 
= 8. The solid line is computed with A fitted with a polynomial of order zero, the dashed line with 
the fit of order one and the dash-dot line with the fit of order two. 

FIG. 4. Modified Lagrange multipliers for the case k{E) = 6 + E/2. The circles indicate values 
computed with eq. (8) , the solid line represents the fit with a polynomial of order zero, the dashed 
line the fit of order one and the dash-dot line the fit of order two. The calculation of the error bars 
is explained in the text. 

FIG. 5. Calculated density of states, relative to the exact density of states, for the case 
k{E) = 6-1- E/2. The solid line is computed with A fitted with a polynomial of order zero, 
the dashed line with the fit of order one and the dash-dot line with the fit of order two. 

FIG. 6. Calculated partition function, relative to the exact partition function, for the case 
k{E) = 6 + E/2. The solid line is computed with A fitted with a polynomial of order zero and the 
dashed line with the fit of order one. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9409001vl 
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This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9409001vl 
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This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9409001vl 
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